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Abstract 

The description of nuclear matter bulk properties in the framework of quark- 
based hadron models in curved unit cells is extended to finite temperature. 
We discuss the analogy and interplay of density and temperature effects in 
this formulation and study, as an application, the temperature and density 
dependence of chiral vacuum properties in the Nambu-Jona-Lasinio (NJL) 
model. 
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The complexity of the strong interactions resides to a large part in their ground state, 
the QCD vacuum, which is the focus of various research directions in nuclear and particle 
physics. One of the currently followed strategies towards uncovering the vacuum structure is 
to study its response to changes in external conditions. Systems with a large baryo-chemical 
potential and/or a high temperature in sufficiently extended volumes of space-time are well 
suited for this purpose and experimentally accessible. 

Indeed, the above conditions are realized in the highly compressed matter of dense stars, 
and they can be produced in relativistic heavy-ion collisions. In both situations the baryon 
density and temperature can become large enough to alter significantly the structure of the 
matter ground state, which is expected to eventually undergo transitions to new phases with 
modified symmetry properties. 

In addition, hadron structure becomes a relevant factor for the matter dynamics at high 
densities. Mainly as a consequence of the changing vacuum properties, it also undergoes 
density and temperature induced modifications, which have been estimated at the mean- 
field level in various hadron models |[T|-§[ . Some of these studies explored a new description 
of nuclear matter bulk properties at zero temperature, which divides the dense environment 
into curved cells of unit baryon number Since high baryon densities, in particular if 

produced in relativistic heavy-ion collisions, are typically accompanied by temperatures of 
the order of the QCD scale, this approach should be extended to finite temperature, and 
this is the purpose of the present letter. 

The crucial feature of the new approach, which generalizes the standard unit cell de- 
scriptions for homogeneous or periodic media, is the form of the unit cell, an intrinsically 
curved three-dimensional hypersphere S^{L) of radius L. The hyperspherical cell provides 
a new setting for the study of finite-density dynamics in hadron models, which is in some 
aspects complementary to more conventional approaches. In addition, it requires much less 
computational effort than traditional unit-cell techniques, but reproduces results of conven- 
tional dense-matter calculations in the Skyrme and Nambu-Jona-Lasinio (NJL) models 
almost quantitatively. 
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The recent generalization of the hypersphere approach to quark-based hadron models 
allows us to perform our study in the standard Nambu-Jona-Lasinio model |^. Our 
focus will be mainly on the combined temperature and density dependence of chiral vacuum 
properties and of the chiral phase transition. The pertinent observables are well studied in 
the hypersphere approach at zero temperature and are expected to be rather strongly 
T dependent, as indicated, for example, by lattice predictions of Tc ~ 140 — 160 MeV for 
the chiral restoration temperature |^. 

In a hyperspherical unit cell S^{L) of radius L, the standard NJL lagrangian takes the 
form 

The g's are the (color triplet and iso-doublet) quark fields with current mass m^, and 
their characteristic four-quark contact interaction is written in terms of the matrices 
r° = I^It, r* = i'j^Ti {i = 1, 2, 3). (The index a G {0, 1, 2, 3} is implicitly summed over in 
the interaction lagrangian above.). In the unit cell context, the geometrical and topological 
properties of 5*^ are mediating interactions with the dense environment by characteristically 
modifying the NJL dynamics: the fiat-space gradients are replaced by their analogues 
Ca (the arrow indicates operation to the left) and a density-induced, isoscalar gauge field 
uJa, the spin connection, emerges. Our notation follows ref. 0, to which we also refer for a 
detailed discussion of the new features of hyperspherical cells in the context of quark-based 
hadron models. 

For the study of the model (P at finite temperature we adopt the standard imaginary 
time formalism 0. The corresponding compactification of the time direction does not 
interfere with the curvature and topology of the hypersphere, which only affects the three 
spatial dimensions. In fact, the descriptions of temperature and density in hyperspherical 
cells do not only coexist naturally but also show direct analogies, which we will discuss 
further below. 

In order to calculate the thermal Gibbs average of any desired Green's function in the 
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given framework, we analytically continue to euclidean times r = ixo and impose the Kubo- 
Martin-Schwinger (KMS) boundary conditions [^]. For the quark propagator, in particular, 
they require Spi^r+P) = — e'^^ Sjs^r), where P is the inverse temperature and ( is the chemical 
potential which sets the average baryon number in the hyperspherical unit cell to one (see 
below). The antiperiodicity in euclidean time is a common property of all (fermionic) finite- 
temperature Green's functions and implies that the model is now defined on a fully compact 
and boundary-less domain, S^{L) x S^{^). 

The NJL quark propagator in Hartree-Fock approximation has a space-time independent 
quark self-energy with both scalar and vector parts, S = —74^^,. Its hypersphere version 
was derived in ref. at T = and can be easily generalized to finite temperature: 

5^3 (x, r) = e^(^^^+^s„.) ^^(^^ ^) ^ ^ ^^^^(^^ ^) _ ^^(^^ ^) ]^ (2) 

The spatial coordinates of the hypersphere are denoted in bold face, i.e. x = {/i, ^,0}, 
f = f{6, 0) is the unit-vector in R^, S is the Dirac spin matrix and the euclidean 7 matrices 
are chosen hermitean, with {'j^,'yu} = '^^iiu- The three amplitudes S'i(/i,r) follow from 
their zero-temperature counterparts by Wick-rotating to euclidean time and imposing 
the KMS condition. In the following we will only need two of them explicitly, namely, 

«.(..^)-^E-We,(A-M E^^jp^ (3) 

and 5*0, which is related to 52 by a time derivative, S'o(/i, r) = m^^ (9^S'2(/i, r). 

These amplitudes represent the propagator in terms of the discrete constituent quark 
spectrum in the unit cell, which is completely determined by the geometrical properties 
of S^{L) and summarized in the generalized momenta /c„ = {2n + 1)/(2L), energies Un = 
\Jk'^ + m? and degeneracies (apart from flavor and color) Dn = 2n{n + 1). 

Analogously, the compactified euclidean time dimension leads to a discretized Matsubara 
frequency spectrum k^^'' of the quarks, which is additionally shifted by the vector part of 
the self-energy and by the chemical potential: 

~kt^ = k^r^ - ^(C + S.) , fcf ^ = (2n + if-. (4) 
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Anticipating the well-known divergencies of quark loops in the NJL model, we also intro- 
duced a soft momentum cutoff A in eq. (^. It has the form of a regularized theta function, 
Ge(A — k) = (^^Q^^i''-^) -)- 1)^^, with e parametrizing the width of the transition region. 

The entire spatial coordinate dependence of the amplitudes 5*0 and S2 is contained in 
the set of functions 

2A;„Lsin(fc„L/i) - cos(fc„L/i)tanf 

SnW = : s„(0) = (5) 

sm/i 

which were introduced in ref. P]. Their index refers to the contribution from the n-th quark 
level. The s„ are orthogonal on and normalized to the corresponding degeneracies: 

— / rf/i sinV Snifi) Smilj) = SnmDn- (6) 
ZTT JO 

With the finite-temperature quark propagator at hand, we can now derive the combined 
temperature and density dependence of the fundamental chiral properties, i.e. of the con- 
stituent quark mass, the chiral order parameter and the pion decay constant. Up to the small 
current masses, the constituent quark mass of the NJL model is dynamically generated as 
the scalar part of the quark self-energy 

E(a;) = S = ig{2TatT [Sp{x,x)Ta] - Ta [Sp{x,x) + S^{x,x)] Ta }. (7) 

(The superscript ± of the propagators distinguishes the coincidence limits from positive and 
negative times.) From eq. we derive the scalar and vector components of the self-energy 
by inserting the propagator (H), 

= -4tg {l + 2n,nf) 82(0) 

2ma °° 00 ^ 

''l + 2n,nj)Y.D^e,{A-K) ^ —-^^ -, (8) 



= 8«(?S'o(0) 



(m) 



= i|f;D„e.(A-w f: . (9) 

n=l m=-oo [kl 'Y + Ui 

(Here, ric and Uf denote the number of quark colors and flavors, respectively.) 
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The Matsubara sums in the above expressions can be rewritten as contour integrals 
and evaluated analytically, with the result 

= ^(1 + 2n,nf) £ ^e,(A - ^{1 - n+iK) - n'{K)}, (10) 
= -f7 E DnQeiA - kn){n-{K) - n+{K)}, (11) 

^ n=l 

where n^{kn) = [e^'^" + are the Fermi occupation number distributions^ with = 

cc;„ ± (C + K)- 

As usual, the dynamical quark mass is obtained as the nontrivial solution of the gap 
equation m{p,T) = rrio + Ss(p, T), which formalizes the Hartree-Fock self-consistency con- 
dition by requiring the value of m in the quark loop of (which appears implicitly in 
the energies Un) to be equal to the dynamically generated constituent quark mass. The 
temperature and density dependence of the solutions to the gap equation will be discussed 
below. 

From the trace of the Hartree-Fock quark propagator we directly obtain the second chiral 
property of interest, namely the scalar quark condensate 

n E 

« qq » = -ztr,,,[5^(0)] = --—f (12) 

which serves as the standard order parameter of the chiral phase transition. (The double 
brackets denote the thermal average of the operator they enclose, and the trace is over 
color and spinor indices.) The quark number density, or equivalently the vector quark 
"condensate", can be similarly expressed as << g^g >>= —12^^,^/ {2g). Since << g^g >> 
is kept fixed in the unit cell at 3p/2 = 3/{2V), the above equation determines the chemical 
potential as a function of the quark mass, the baryon density and the temperature, i.e. 
C = C{m,p,T). 



^Note that the zero-temperature limit of eqs. ( p!o| , p!l| ) leads back to the expressions of ref. 
(with three valence quarks in the cell). The subsequent L ^ 00 limit reproduces the familiar 
zero-density results of the NJL model Q, as shown in ref. (P). 



The remaining chiral quantity of interest is the pion decay constant /tt- We obtain it 
from the Bethe-Salpeter amphtude X°'{^jy'iP4^) of the pion in its rest frame {i.e. at p = 0), 
which is proportional to the one-pion matrix element of the axial current and thus to /yr: 



« 0\Al{x)\n\p,) » =-tT,,cj [7m75^x'(x,x;p4)] = -pJJ'"'^^^=. (13) 



2ujpV ^2ujpV 

Due to the periodicity in the euclidean time direction, the pion energy (which is later analyti- 
cally continued back to the pion mass) becomes a bosonic Matsubara frequency P4 •* = 2l7r/(3. 
The trace is over color, flavor and spinor indices and connects the external quark lines of 
the Bethe-Salpeter amplitude to the axial vector vertex. 

The solution of the NJL Bethe-Salpeter equation for the pion in ladder approximation 
is, apart from the external quark propagators, just the pseudoscalar pion-quark Yukawa 
vertex, with the coupling determined by the Goldberger-Treiman relation at the quark level, 
g-TTqq = TTT'/ fir- Adapted to finite temperature it reads 

777 r 

X''{x,y-p,) = -- / dfi{z)e-'P^^^Sp{x,z)T^j,Sp{z,y) +0{pl). (14) 

Jtt 

The integration range in eq. ([l^) is S^{L) X S\^) with the measure d^{z) = 
Li^ siv? iismO dn dO dcj) dTz- From eq. (1T3| ) we now obtain 

fl = AtTrn.nf— ^ dfi sinV dre-'^^^ 8^ {S2ii^,T) 82(1^, -r)}, (15) 

P4 Jo Jo 

and after a partial integration in the r integral (the surface term vanishes, due to the 
periodicity properties of both 5*2 and the exponential) we insert the explicit form (^) of 
the propagator amphtude 5*2. The integration over the geodesic hypersphere coordinate fi 
is immediately performed by using the orthogonality relation (|^), and the integral over r 
combines the two independent Matsubara sums into one, 

1-277 rl3 



/ (i/i sin fi / dre ^ S'2(/U, r) S'2(/U, — r) 
JO Jo 



n m'^ °° °° 1 1 
'^^E^n0e(A-A;„) E m 17^ m ' (16) 
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which can again be evaluated exphcitly in terms of contour integrals. After continuing back 
to Minkowski energies and specializing to the chiral limit, we arrive at our final expression 
for the square of the pion decay constant: 

= E ^ 0e(A - K) {l - n^iK) - n-{K)} . (17) 

Note that the explicit temperature dependence (the factor in curly brackets) both of this 
expression and of the scalar quark self-energy (|10|) are identical. The overall T dependence 
differs, however, due to the implicit T dependence of the u;„. Nevertheless, the factor of 
ensures the decoupling of the pion at the chiral phase transition. 

The final expressions for the gap equation and the chiral vacuum properties have the sim- 
ple and intuitive form which is typical for hypersphere results, despite the perhaps somewhat 
unfamiliar formalism used in their derivation. In fact, the density-induced modifications of 
the fiat-space results affect primarily the constituent quark spectrum and not the general 
form of the expressions. This makes the transition to finite temperature particularly straight- 
forward in the NJL model, where a constant mean field (the constituent quark mass) breaks 
chiral symmetry. It can be "heated" as in standard perturbation theory, since it is explicitly 
generated by loops. This is in contrast to the Skyrme model, for example, where thermal 
fluctuations around the spatially varying soliton field would have to be taken into account. 

To prepare the quantitative discussion of our results, we fix the model parameters at the 
same values as in ref. 0: ^ = 4.08 GeV-^ A = 700 MeV and e = 2.4. The behavior of 
the chiral order parameter as a function of the temperature is shown in Fig. 1 for different 
values of the density. At small T the condensate is only weakly temperature dependent, 
even at rather high densities. With increasing temperature it starts to drop rather steeply 
on the approach to the the critical temperature of the chiral restoration transition. As 
one would expect, Tc itself decreases monotonically with density. 

Figure 2 displays the temperature dependence of the pion decay constant, again with 
the density as a parameter. As expected, f-,, reaches zero at the same (density dependent) 
critical temperature as the quark condensate, confirming that the pion decouples from the 
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vacuum when chiral symmetry is restored. However, f.^ starts to decrease significantly at 
somewhat lower temperatures than the quark condensate, exactly as was observed in the 
conventional chemical potential approach [Q]. 

Since the description of density and temperature in hyperspherical cells is built on the 
compactification of space and time directions, respectively, the explicit density and tem- 
perature dependence of the chiral observables factorizes. The implicit coupling through 
the solution m(p, T) of the gap equation leads, nevertheless, to a non-trivial interplay be- 
tween temperature and density. It is therefore remarkable that the temperature and density 
dependence of the observables follows rather closely the one obtained in the conventional 
treatment 0, even though it can be traced to specific features of the hyperspherical cell. 
This seems to be a recurrent theme in hypersphere calculations, already encountered in the 
studies of the Skyrme and NJL models at zero temperature. 

We conclude with several comments. Our study of the combined temperature and density 
dependence of chiral properties in the Nambu-Jona-Lasinio model is the first application of 
the hypersphere approach at finite temperature. The extension to finite temperature fits 
naturally into the hypersphere framework and extends its range of applicability a step farther 
towards (if only qualitative) contact with experiment. 

Our results resemble rather closely those of the conventional chemical potential approach 
to the NJL model, which approximates hot and dense matter as a constituent quark gas. 
Some thermodynamical quantities are sensitive to this neglect of the quark clustering into 
color singlets [|ll|, and the hypersphere description does not resolve this problem, despite 



its inclusion of, on average, three quarks into a unit cell. The latter does not lead to 
confinement or clustering of the quarks, since the periodicity of the hypersphere coordinates 
allows them to enter neighboring cells by travelling more than once around the sphere, in 
complete analogy with the effect of periodic boundary conditions in standard unit cells. 

Our treatment of the temperature dependence in the unit cell is in some sense minimal, 
since one could have allowed a temperature dependence of the unit-cell metric in addition 
to the thermodynamics of the quarks. It is conceivable, for example, that an extended 
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variational approach would find a different cell shape favorable at high temperatures, leading 
to an analog of a Landau transition in crystals. This is not likely, however, since the cell 
metric has its origin mainly in the chiral symmetry properties of the model dynamics, which 
do not change at finite temperature. 

Finally, it is interesting to note that the realization of finite density and temperature 
are remarkably similar in hyperspherical unit cells. In particular, both are described by 
compactifying space-time directions, leading to discrete momenta and energies, respectively, 
with their scale set by p and T. (The time direction stays flat, of course, since a one- 
dimensional circle has no intrinsic curvature.) 

Since there are many similarities in the effects of finite baryon density and temperature, 
for example on phase transitions and on hadron properties, it is tempting to speculate that 
a closely analogous description of both would make their common aspects more transparent 
and that the hypersphere approach gives a hint in this direction. 

This work was supported by the U.S. Department of Energy under grants no. DE-FG02- 
93ER-40762 and DE-FG02-88ER-40388. 



11 



REFERENCES 

[1] E. Shuryak, The QCD Vacuum, Hadrons and the Superdense Matter (World Scientific 
Pub. Co., Singapore, 1988) and references therein 

[2] I. Klebanov, Nucl. Pliys. B 262 (1985) 365; L. Castillejo, P.S.J. Jones, A.D. Jackson, 
J.J.M. Verbaarschot and A. Jackson, Nucl. Phys. A 501 (1989) 801 

[3] T. Hatsuda and T. Kunihiro, Phys. Rev. Lett. 55 (1985) 158;V. Bernard, U.-G. Meissner 
and I. Zahed, Phys. Rev. D 36 (1987) 819; T. Hatsuda and T. Kunihiro, Phys. Rep., 

to be published, and references therein 

[4] N.S. Manton and P.J. Ruback, Phys. Lett. 181 B (1986) 137; N.S. Manton, Commun. 
Math. Phys. Ill (1987) 469; A.D. Jackson, NORDITA preprint 87-38 N (1987) 

[5] A.D. Jackson, A. Wirzba and L.Castillejo, Phys. Lett. 198 B (1987) 315, Nucl. Phys. 
A486 (1988) 634; H. Forkel, A.D. Jackson, M. Rho, C. Weiss, A. Wirzba and H. Bang, 
Nucl. Phys. A 504 (1989) 818; A. Wirzba and H. Bang, Nucl. Phys. A 515 (1990) 571 
and references therein 

[6] H. Forkel, Phys. Lett. 280 B (1992) 5; University of Maryland preprint #94-037 (1994), 
to appear in Nucl. Phys. A 

[7] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 124 (1961) 246, 255 

[8] J.B. Kogut, D.K. Sinclair and K.C. Wang, Phys. Lett. 263 B (1991) 101 

[9] J.I. Kapusta, Finite-Temperature Field Theory, Cambridge University Press, Cambridge 
(1989) 

[10] R. Kubo, J. Math. Soc. Japan 12 (1957) 570; P.C. Martin and J. Schwinger, Phys. Rev. 
115 (1959) 1342; R. Haag, N.M. Hugenholz and M. Winnink, Commun. Math. Phys. 5 

(1967) 215 

[11] K. Kusaka, Phys. Lett. 269 B (1991) 17 



12 



FIG. 1. — << qq >>^/^ as a function of temperature for five values of the baryon 
density, p/pQ = {1.58,1.34,0.89,0.46,0.01}, in units of the nuclear matter saturation density 
po = 0.17 fm~^. (The uppermost curve corresponds to the lowest density, and the ones below 
are arranged in order of increasing density.) 

FIG. 2. The pion decay constant as a function of temperature, parametrized by the same five 
values of the density as in Fig. 1. 
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